Method for detecting the need for calibrating probes

ABSTRACT

An in situ probe failure detection method for probes which deliver an electric signal on account of a change of the concentration of a substance to be measured. An example of such a probe failure may be caused by, for example, fouling of the membrane as the relative change in signal may be interpreted as a change in the concentration of the substance to be measured. More particularly, the method involves detecting the need for calibrating probes, which probes deliver an electric signal on account of a change of the concentration of a substance to be measured, characterized by: (a) subjecting the probe to a stepwise change of the concentration of the substance to be measured; (b) calculating the value of the probe response time constant t c  of the obtained measurement curve with the help of the following main equation (1): t c  dy/dt=-y+α 0  +α 1  t wherein t=time (time), t c  =probe response time constant (time), y=substrate probe signal (mass volume -1 ), α 0  =substrate concentration at start of response (mass volume -1 ), α 1  =rate of change of substrate concentration (mass volume -1  time -1 ), and (c) comparing the calculated t c  values of at least two different measurements with each other, and in case the calculated t c  values do not match, calibrating the probe.

The invention relates to a method for detecting the need for calibrating probes which probes deliver an electric signal on account of a change of the concentration of a substance to be measured.

BACKGROUND OF THE INVENTION

A general problem with probes as indicated above, for instance, dissolved oxygen measuring probes (DO probes) is that probe failure caused by, for example, fouling of the membrane may be interpreted as a change in the concentration of the substance to be measured.

Normally this problem is solved by calibrating the probe on a regular time basis. However this solution is not without problems, because by periodically calibrating the probe, the calibration can be either too late or too early. When the calibration is done too early it costs too much working-hours and if the calibration is too late the probe had produced false values during an unknown period.

SUMMARY OF THE INVENTION

The above-indicated problem has been solved with the invention presented below. It has been found that the probe response time constant can be calculated periodically and that this parameter can be used for probe diagnosis.

The method is based on the fact that the response of the probe after an instant change in concentration can be modelled as a first order model combined with a ramp function: ##EQU1## wherein t=time (time), t_(c) =probe response time constant (time), y=substrate probe signal (mass volume⁻¹), α₀ =substrate concentration at start of response (mass volume⁻¹), α₁ =rate of change of substrate concentration (mass volume⁻¹ time⁻¹). The measured response of the probe is used to calculate the constants in this equation. The probe response time constant (t_(c)) from this model can be used for probe diagnosis. As long as the probe response time does not change in consecutive measurements of the probe response after an instant change in concentration the probe is accurate. As soon as the probe response time is changed during consecutive measurements the probe should be recalibrated. In the method according to the invention the concentration of the substance to be measured is changed periodically and each time after such a change the probe response time constant is calculated. As long as the probe response time constant doesn't vary there is no problem. As soon however as the probe response time constant changes there is a probe failure and the probe has to be cleaned and calibrated. By using the method according to the invention described a new calibration is always done in time.

To prove the method, Applicant has carried out research with a continuous flow-through respiration meter (EP-B- 257057). The continuous flow-through meter measures the dissolved oxygen (DO) concentration at the inlet and at the outlet of a closed respiration chamber through which the sludge is pumped continuously Spanjets H. and Klapwijk A., (1990). On-line meter for respiration rate and short-term biochemical oxygen demand in the control of the activated sludge process. Advances in Water Pollution Control, Proceedings of the 5th IAWPRC Workshop held in Yokohama and Kyoto, Japan, 26 Jul.-3 Aug. 1990!. The respiration rate is calculated from the difference of the two DO measurements.

The electrochemical oxygen sensor used in the continuous flow-through respiration meter is based upon the electrochemical reduction of oxygen in an amperemetric cell. The electrodes are immersed in an electrolytic solution which is separated from the bulk solution by a semi-permeable membrane. The electrode signal is determined by the diffusion of dissolved oxygen (DO) from the bulk through the membrane. The probe-meter combination will be referred to as "DO-probe".

The respiration meter as described by Spanjets and Klapwijk (1987): Measurement of the toxicity of KCN and some organic compounds for the activated sludge using the Wazu-respiration meter; In: Kolaczkowski S. T. and Crittenden B. D. (eds.) Management of hazaradous and toxic wastes in the process industries; Elsevier, and Spanjars and Klapwijk (1990) loc.cit., has proven its usefulness and its reliability. The characteristic of this meter is that the DO concentration in the inlet and outlet is measured with one single probe, fixed at one opening of the chamber, by alternating the flow direction through the chamber (FIG. 1).

When the DO probe is replaced by another probe or sensor, such as a nitrate or a chlorine probe, the method according to the invention is also applicable to this kind of probes or sensors.

An aspect of this principle is that, for instance, a steady state DO reading is only available when the DO probe signal, after changing the flow direction, has reached its end value. Hence the maximum measuring frequency for the DO concentration is limited by the response time of the DO probe. Till now, in practice, the measurement of the DO takes place when the signal has reached about 95% of its ultimate value. The time needed to attain this 95% response (response measuring period) is based on experience. As both the inlet- and outlet concentration may be off the real value, a systematic error of about 10%, may occur in the calculated respiration rate. Furthermore s dirty or fouled membrane and consequently s slower DO response may result in even s higher error contribution.

As elucidated below s method is proposed to estimate the steady state DO concentration from the probe response curve, even though the full response has not yet been attained.

Furthermore the first order probe response time constant will be estimated after each response, which can be used for probe performance diagnosis. The method according to the invention is based on a least squares fit to a first order probe response model. Simulations and batch experiments are used to demonstrate the effectiveness of the method.

MODEL OF THE DO-PROBE RESPONSE

As indicated above in is assumed that the DO probe signal can be modelled by a first order dynamic system. Additionally, it is assumed that the change of the real DO concentration during the response measuring period will change along a linear slope. The response of the DO probe after a change in the flow direction can then be considered as a combined step and ramp response of a first-order system (FIG. 2).

In accordance with the invention the probe signal, for instance the DO probe signal y is modelled by a first order dynamic system: ##EQU2##

Where the input c is the real DO concentration which is assumed to be a ramp function with offset α₀ and slope α₁ :

    c=α.sub.0 +α.sub.1 t                           (2)

Combining (1) and (2) gives the response function: ##EQU3##

There is an analytical solution of (3) that can be written in the form: ##EQU4## where: y(0)=the initial condition of the probe signal at time t=0.

    u=(α.sub.0 +α.sub.1 t)/t.sub.c

Inserting the value of u gives the DO response as a function of time:

    y(t)= (y(0)-α.sub.0 +α.sub.1 t.sub.c)e.sup.-t/tc !+ α.sub.0 +α.sub.1 t!- α.sub.1 t.sub.c !(5)

It is recognized that the first term of (5) converges to 0 after some time, determined by the probe response time. The last term of (5) is then seen to show the lag between the real DO concentration and the probe response after the transient. FIG. 3 shows the result of a simulation using (5) for different realistic values of the parameters. The simulation demonstrates that, when the real DO concentration changes along a linear slope, the ultimate value of the DO signal (y) always deviates α₁ t_(c) from the real DO concentration (c).

ESTIMATION OF DO CONCENTRATION AND PROBE RESPONSE TIME

All the three parameters t_(c), α₀ and α₁ in equation (3) are unknown and have to be estimated from the DO probe step response. Since the DO measurements are obtained in time discrete form it is natural to approximate the time derivative in (3) with finite differences. Here Applicant has chosen the trapezoidal rule for integration: ##EQU5## where: k=1, 2, . . . , K; K typically 10 or 15 h=sampling interval DO concentration, typically 1 or 2 s.

The parameter α₁ can be estimated separately. Here Applicant has considered the last DO measurement of two consecutive step responses, called y_(K) (m-1) and y_(K) (m). The DO concentration slope is then estimated from: ##EQU6##

Now there remain two parameters, t_(c) and α₀, to be estimated. If equation (1) is applied to all K measurements (but the first) in one step response the following equations are obtained:

    t.sub.c (y.sub.2 -y.sub.1)-hα.sub.0 =1/2h(-y.sub.2 -y.sub.1 +3hα.sub.1)

    t.sub.c (y.sub.3 -y.sub.2)-hα.sub.0 =1/2h(-y.sub.3 -y.sub.2 +5hα.sub.1)

    t.sub.c (y.sub.k -y.sub.k-1)-hα.sub.0 =1/2h(-y.sub.k -y.sub.k-1) +(2k-1)α.sub.1)                                     (8)

This can be written in a vector form: ##EQU7## or written in a more compact form:

    φθ=y                                             (10)

The unknown parameter vector θ can be estimated by the least-squares method (Astrom K. J. and Wittenmark B., (1984) Computer Controlled Systems, Theory and Design, Prencetise-Hall, Inc., New York):

    θ=(φ.sup.T φ).sup.-1 φ.sup.T y           (11)

The estimate θ makes it possible to calculate the DO concentration at the end of the sampling interval:

    c(m)=α.sub.0 +α.sub.1 Kh                       (12)

This is used for a new estimate of α₁ : ##EQU8## c(m-1) is the calculated real DO from the foregoing response. Then a new vector θ is calculated from the observations y. This procedure is repeated until a preset error criterion is met.

This method was tested by means of two simulations and an experiment. The application of the method was demonstrated in another experiment where the probe membrane was fouled intentionally.

MATERIALS AND METHODS

In the simulations, executed with SIMON Elmqvist H. et al., (1986) Simnon, User's guide for MS-DOS computers, SSPA Systems, Sweden!, measuring data were generated with the use of the first order model (5) and chosen values of c and t_(c). Then the method proposed was used to recalculate c and t_(c). Any deviation of the result from the original c and t_(c) emanates from the assumption that the time derivative of (3) can be approximated by the trapezoidal integration method.

In the laboratory experiments a respiration meter RA1000, Manorharm, equipped with a WTW dissolved oxygen meter (modified model OXY-219/R with sensor model E090), was connected to an aerator, the total system having a content of 1.5 liter of activated sludge. The sludge was sampled from a nitrifying activated sludge plant. The respiration meter which works according to the principle described by Spanjars and Klapwijk (1990 loc.cit.) was operated in a mode such that the alternating DO signal could be sampled. Temperature and pH were kept constant at 20° C. and 7.5 respectively.

In Example 1 the assumption of the first order model was tested. Therefore a known amount of ammonium was added to the sludge and the DO signal was recorded. The DO concentration at the inlet and at the outlet were either calculated from 2-4 values at the end of the response measuring period (averaging method) or calculated according to the method of the invention (estimation method). From the DO concentrations, the respiration rate and the mass of oxygen used by the nitrifiers for the oxidation of the unit of mass nitrogen (O/N ratio) were calculated respectively. The O/N ratio was compared for different situations.

In Example 2 the procedure to detect probe failure from an increasing t_(c) was tested. Probe fouling was imitated by covering a part of the probe membrane with ball-bearing grease. Ammonium was added to the sludge before and after contamination the membrane. The effect on t_(c) and on the O/N ratio was studied.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1: Scheme of the one-probe continuous respiration meter.

FIG. 2: Combined step and ramp response of the DO probe after changing the direction of the sludge flow. First-order linear system assumed.

FIG. 3: Simulation, using equation (5), of the measured DO (y) after reversing the sludge flow direction i.e. after a step change in the DO (c) followed by a time linear change of c. Probe time response constant t_(c) =0.10 min.

FIG. 4: Simulation, using (5), of the measured DO concentration (points) during the operation of the respiration meter in comparison to the inlet and outlet-DO (solid lines). h=2 s, t_(c) =0.08 min, t_(r) =0.5 min, α₀ 32 4 g m⁻³ and 2 g m⁻³ for respectively the inlet and the outlet, α₁ =0.4 g m⁻³ min⁻¹ and 0.5 g m⁻³ min⁻¹ for respectively the inlet and the outlet. Standard deviation noise: 0.008 g m⁻³.

FIG. 5: Result of the estimation method applied on the simulated DO probe values (FIG. 4). Calculated values (points) compared with originally simulated values (solid lines). a. Probe response time constant (t_(c)). b. Real DO concentration (c).

FIG. 6: Simulation of a batch experiment. Addition of 0.010 g ammonium-nitrogen to 1.5 liter sludge. Monod constants: K =1.5 g ammonium-nitrogen per m³, maximum respiration rate=65 g m⁻³ h⁻¹ ; Probe response time constant (t_(c)): 0.1 min; standard deviation noise: 0.01 g O₂ m⁻³ ; Response measuring period (t_(r)): 30 s.

FIG. 7: Result of the estimation method for c and t_(c) applied on the simulated DO probe values (FIG. 6). Estimated values (points) compared with the originally simulated values (solid lines). a: Probe response time constant (t_(c)). b. Real DO concentration (c).

FIG. 8: DO-probe signal in a batch experiment: addition of 5 mg N as ammonium to activated sludge. a: full experiment; b: first 1.5 minutes of the experiment.

FIG. 9: t_(c) calculated from the DO-probe measurement. a: without correction for dead time; b: with correction for dead time.

FIG. 10: Respiration rate calculated from the DO probe signal in FIG. 8. Comparison of the averaging method and the estimation method according to the invention.

FIG. 11: Effect of the contamination of the probe membrane. a: Probe response time constant (t_(c)). b: Respiration rate; comparison of the averaging method and the estimation method according to the invention. c: Part of the DO probe signal (points) before and after contamination (respectively left and right), in comparison with the estimated DO.

FIG. 12: Addition of 0.005 g N (as ammonium) to 1.5 l of activated sludge, before and after contamination of the probe membrane. a: Respiration rate calculated from the DO concentration. Comparison of the averaging method and the method according to the invention. b: Probe response time constant (t_(c)).

EXAMPLE 1 Verification of the method Test 1: Linear change of the real DO concentration (simulation)

For given constant values of α₁, h and t_(c), and an initial value of α₀, equation (5) was used to simulate the course of the measured DO during a short period of respiration measurement (FIG. 4). In this simulation it was assumed that the respiration rate was decreasing linearly. Normally distributed noise was added to the calculated value of the DO.

The simulated DO values were used to estimate t_(c) and c employing the equations (7) through (13). FIG. 5 shows the result.

In case of an optimal estimation procedure, the resulting t_(c) and c would have been the original points. The figure shows that t_(c) (mean value 0.079±0.002 min) and c are estimated well.

Test 2: Non-linear change of the real DO concentration (simulation)

A batch experiment with the continuous flow-through respiration meter was simulated to test the assumption that, within the response measuring period, the change of the real DO concentration can be accurately represented by a linear relationship. In this experiment an amount of ammonium was added to the aerator connected to the respiration meter. The simulation model included Monod kinetics and a set of dynamic mass balances on DO and substrate for the aerator and the respiration chamber. The parameters were chosen from a fit of the simulation model on the experimental results of Spanjars and Klapwijk (1990 loc.cit.). In this test, unlike test 1, the change of the simulated DO is not linear any longer. Even so, in the calculation of t_(c) and c, it is assumed to be linear within the response measuring interval.

FIG. 6 shows the simulated DO concentration at the inlet and at the outlet of the respiration chamber as well as the simulated response of the DO probe.

The simulated DO probe measurement values were used to estimate t_(c) and c with the equations (9) through (13). FIG. 7 shows the result.

As in test 1, t_(c) and c are estimated well (FIG. 7). t_(c) does not depend on the magnitude of the response. The estimated DO is very close to the originally simulated DO.

From the results of test 1 and 2 it is concluded that the method proposed here allows the estimation of the first order response time constant and the real DO concentration.

Test 3: Experimental verification of the first order model

To test the assumption of the first order model, the experiment described and simulated in test 2 was carried out in the laboratory. In three experimental runs, the response measuring period (twice 15 s and once 20 s) was too short to attain the steady state response. Therefore this experiment was well suited to demonstrate the capability of the method to estimate the real DO concentration in the inlet and in the outlet of the respiration chamber. FIG. 8 shows the DO-probe signal.

FIG. 9a shows the estimated t_(c) for both the upward and the downward response. t_(c) from the upward response is higher than t_(c) from the downward response. Apparently, the model applied is not completely valid for both the upward and the downward response or both. After careful examination of the shape of the response curves (FIG. 8b) in relation to the geometry of the respiration meter (there was a dead space in one flow direction) it was concluded that a dead time (t_(d)) had to be incorporated in the model (3): ##EQU9## t_(d) depends on the flow direction. In the experiments reported here, the values were found by trial and error to be 2 s and 1 s for the measurement of the inlet DO concentration and the outlet DO concentration respectively. FIG. 9b shows that when the dead time is accounted for, both responses produce almost similar t_(c) 's which also show less variation. As in test 2, t_(c) is independent of the response amplitude.

The DO concentration at the inlet and at the outlet of the respiration chamber were calculated by averaging three measurements at the end of each response (averaging method) and by the method according to the invention. Next, the respiration rate was calculated as described by Spanjers and Klapwijk (1990 loc.cit.). The result is shown in FIG. 10.

The method according to the invention results in higher values for the respiration Pate, especially at a higher rate when the difference between inlet and outlet DO concentration is greater. It is difficult to verify the correctness of the absolute value of the respiration rate. Therefore, from the respiration rate, the total amount of oxygen additionally used to the endogenous oxygen consumption, defined as Short Term BOD (BOD_(ST)), was evaluated for both results. The BOD_(ST) was calculated from the area under the respiration curves. This value divided by the amount of nitrogen added represents the mass of oxygen used by the nitrifiers for the oxidation of one unit of mass nitrogen (O/N ratio). Table 1 summarizes the results for the same experiment repeated with different response measuring periods.

                  TABLE 1     ______________________________________     Comparison of the average method (average of three     measurements at the end of the response measuring period)     and the method according to the invention for different     response measuring periods (t.sub.r). Theoretical value O/N:     4.57; expected value 4.4.                O/N          method     h    t.sub.r                averaging    according to     (s)  (s)   method       the invention                                      reference     ______________________________________     1    15    3.91           4.26     this investigation     1    15    3.74           4.01     this investigation     2    20    4.29           4.39     this investigation     2    30    4.02           4.04     this investigation     2    30    4.33           4.38     this investigation     --   30    4.36           --       Spanjers and Klap-                                        wijk, 1990 loc.cit.     --   --    --        4.33 --       Sharma and Ahlert,                                        1977*     ______________________________________      *Sharma B. and Ahlert R.C. (1977), Nitrification and nitrogen removal.      Wat. Res. 11, pp. 897-925

From table 1 two conclusions can be drawn.

Firstly, for a short response measuring period, the averaging method yields lower values for the O/N ratio while the estimation method according to the invention produces equal values compared to a large response measuring period. This means that the estimation method produces a better estimate of the real DO concentration. Two experiments yield a low O/N value from the estimation method compared to the value in the literature.

Secondly, the negligible difference between the two methods at a response measuring period of 30 s indicates that the effect of a changing real DO concentration on the probe response is marginal. For, in the determination of the O/N ratio there is a sharp change of the DO at only two different occasions. Another reason for the small difference is that the shift in both the endogenous as well as in the maximum respiration rate reduces the impact on the total amount of oxygen calculated from the area under the curve. Nevertheless, in kinetic experiments, when transients in the respiration rate and thus in DO concentration become important Ossenbruggen P. et al., (1991) Designing experiments for model identification of the nitrification process; Watermatex 1991!, the effect of the changing real DO will become stronger in the probe response.

From test 3 it is concluded that, for incomplete responses, the method according to the invention produces a better estimate of the real DO concentration and consequently a better result for the calculated respiration rate. Because of the actual geometry of the respiration meter, a correction in the response model for dead time is required.

EXAMPLE 2 Detection of probe fouling using the probe time constant.

In batch experiments is was attempted to detect probe failure from the estimated t_(c). Therefore, the probe membrane was contaminated intentionally by covering part of the membrane with a thin layer of grease (approximately 10-25% of the surface).

First it was established that the contact of approximately the same amount of grease with the sludge in the aerator had no impact on either the DO concentration or the respiration rate. This means that the substance is not biodegradable or toxic on a short term and that it should neither provoke a change in the respiration rate when it is applied to the membrane surface.

Two experiments were carried out. In the first experiment an amount of grease was applied twice to the membrane, while the respiration meter was measuring the endogenous respiration rate. In FIG. 11a the estimated t_(c) is shown and the applications of grease are indicated. The figure shows that a greater amount of grease results in an increasing t_(c) meaning that the probe becomes slower. Consequently, the DO concentration calculated from the end values of the response measuring period (averaging method) is erroneous and so is the calculated respiration rate (FIG. 11b). The DO from the method according to the invention and so the calculated respiration rate is affected to a lesser extent by the deterioration of the probe: after the first treatment the rate decreases from 20 to 15 gm⁻³ h⁻¹ ; after the second treatment it remains on this level, although the noise has increased. This decrease may partly be caused by another incident, for example the real endogenous respiration rate may have decreased during the time that the membrane was being contaminated. This is supported by the check at the end of the experiment where the membrane was replaced by a new one and the respiration rate was found to be 18 mg l⁻¹ h⁻¹.

FIG. 11c illustrates the effect of the contamination on the individual response curves. In this figure the estimated DO according to the invention is also indicated.

In the second experiment ammonium was added twice to activated sludge in the endogenous phase of respiration. The first addition is already reported in the previous section (FIG. 10 and table 1). The second addition was done after contamination of the membrane, FIG. 12a shows the respiration rate calculated from the DO concentration, whereas t_(c) is plotted in FIG. 12b. The O/N ratio was calculated from the respiration rate (table 2).

                  TABLE 2     ______________________________________     Effect of the contamination of the probe membrane on     the calculated O/N value. Comparison of the averaging     method (averaging three measurements at the end of the     response measuring period) and the method according to     the invention. h = 1 s, t.sub.r = 15 s, theoretical O/N     ratio: 4.57.                 O/N                            method according                 averaging method                            to the invention     ______________________________________     before contamination                   3.91         4.26     after contamination                   3.04         4.28     ______________________________________

From table 2 it can be concluded that the O/N ratio evaluated from the method according to the invention is not affected by the probe deterioration while the O/N ratio from the averaging method (already lower because of the short response measuring period, see table 1) is decreased.

CONCLUSION

As indicated above the purpose of the investigation of Applicant was to develop a method to detect the need for calibrating a probe and to improve the measurement of substances with probes which deliver an electric signal such as DO-probes. In the investigation a respiration meter was used with one and the same DO-probe at the inlet and at the outlet of a respiration chamber. Inherent in the principle of the respiration meter is that the DO-probe is repeatedly subjected to step changes in the DO concentration. This fact is employed, for the improvement of the DO measurement, by fitting a first order response model to each measured response, which provides an estimation of the real DO concentration and of the first order probe response time constant. The change of the real DO concentration during the response measurement is accounted for in the model.

The simulations show that the method allows the calculation of c and t_(c). On the condition that the probe signal can be modelled by a first order dynamic system, a correct estimation of c and t_(c) can be obtained from experimental data. The type of the response is chiefly determined by the manufacturer of the DO meter. However the response is also partly determined by conditions of the respiration meter in which the probe is mounted. The experiments show that, in this case, the assumption of a first order response is reasonable, provided a dead time is accounted for. This dead time is fixed by the respiration meter and can be determined exactly by measuring the dead space. The experiments also show that it is reasonable to assume that the change of the real DO concentration within one response measuring period can be approximated by a linear relationship. This can be clearly illustrated by looking at an enlarged portion of FIG. 8a.

The DO concentration c can be estimated from only a part of the probe response, where the signal has not yet reached its steady state. The advantage, in the respiration measurement, is that the measuring frequency of the DO's at the inlet and at the outlet and herewith the measuring frequency of the respiration rate can be increased.

t_(c) indicates malfunction of the probe even in case of severe deterioration (FIG. 11a), when the probe becomes very slow and the estimation of c likely becomes unreliable.

The method according to the invention, for the estimation of t_(c) and c, can also be used, in combination with the measuring technique, in the measurement of the DO concentration in an activated sludge reactor. On the condition that the respiration meter is installed close to the reactor, the DO concentration in the inflow of the respiration meter is equal to the concentration in the reactor at the sampling point Kim C. W. et al, (1991), Continuous determination of oxygen transfer coefficient with an on-line respiration meter. Proceedings: 3rd IAWPRC regional conference on development and water pollution control, November 20-24, Shanghai!. Besides the measurement of the DO concentration and the optional measurement of the respiration rate this technique, through the estimated t_(c), provides a continuous diagnosis of the probe condition.

As appears from the Examples reported above, the method according to the invention provides a reliable estimate of the real DO concentration from the probe response signal when the probe is subjected, in the respiration measurement, to a repeated, stepwise changing DO concentration. As a result, the reliability of the calculated respiration rate is improved.

The first order probe response constant t_(c), evaluated from the probe signal is 8 useful indicator for detecting fouling of the probe membrane. 

I claim:
 1. Method for detecting the need for calibrating a probe, which probe delivers an electric signal on account of a change of the concentration of a substance to be measured, comprising:(a) subjecting the probe to consecutive changes of the concentration of the substance to be measured by providing the probe alternately with at least two different concentrations; (b) calculating the value of the probe response time constant t_(c), the constants α₀ and α₁ of a probe response curve after the instant change of the concentration with the following main equation: ##EQU10## wherein t=time (time) t_(c) =probe response time constant (time) y=substance probe signal (mass volume⁻¹) α.sub. = substance concentration at start of response (mass volume⁻¹) α₁ =rate of change of substance concentration (mass volume⁻¹ time⁻¹) and (c) comparing the calculated t_(c) -values of at least two different measurements with each other and--in case the calculated t_(c) -values do not match--calibrating the probe.
 2. Method according to claim 1, wherein the probe is a dissolved oxygen (DO) measuring probe.
 3. Method according to claim 2, wherein the consecutive changes of the DO concentration is carried out by providing the DO measuring probe alternately with at least two flows having a different DO concentration.
 4. Method according to claim 2, wherein the consecutive changes of the DO concentration is carried out by alternately reversing the flow direction of a reacting DO containing medium, flowing through a closed chamber provided with one single probe at only one opening of said chamber.
 5. Method according to claim 4, wherein the DO containing medium is activated sludge.
 6. Method for calculating the real value of a probe, which probe delivers an electric signal on account of a change of the concentration of a substance to be measured, comprising:(a) subjecting the probe to consecutive changes of the concentration of the substance to be measured by providing the probe alternately with at least two different concentrations; (b) calculating the real value of the probe which is equal to α₀ +α₁ t of a probe response curve after the instant change of the concentration with following main equation: ##EQU11## wherein t=time (time) t_(c) =probe response time constant (time) y=substance probe signal (mass volume⁻¹) α.sub. = substance concentration at start of response (mass volume⁻¹) α₁ =rate of change of substance concentration (mass volume⁻¹ time⁻¹). 